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Abstract 

We consider the transverse-momentum (g-r) distribution of a diphoton pair produced 
in hadron collisions. At small values of qt , we resum the logarithmically-enhanced 
perturbative QCD contributions up to next-to-next-to-leading logarithmic accuracy. 
At intermediate and large values of qt , we consistently combine resummation with the 
known next-to-leading order perturbative result. All perturbative terms up to order 
«g are included in our computation which, after integration over qt , reproduces the 
known next-to-next-to-leading order result for the diphoton pair production total cross 
section. We present a comparison with LHC data and an estimate of the perturbative 
accuracy of the theoretical calculation by performing the corresponding variation of 
scales. In general we observe that the effect of the resummation is not only to recover 
the predictivity of the calculation at small transverse momentum, but also to improve 
substantially the agreement with the experimental data. 
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1 Introduction 


Diphoton production at hadronic colliders is a very relevant process, both from the point of view 
of testing the Standard Model (SM) predictions [H El [3l HJ |5l | 6 ] as for new physics searches. 

While a general understanding of QCD processes at hadronic collisions poses a serious chal¬ 
lenge due the complicated environment, direct or prompt photons provide an ideal probe since they 
constitute a theoretically and experimentally clean hnal state. From the theory side, because they 
do not have QCD interactions with other hnal state particles and, from the experimental side, be¬ 
cause their energies and momenta can be measured with high precision in modern electromagnetic 
calorimeters. 

Besides purely QCD-related considerations, diphoton hnal states have played a crucial role in 
the recent discovery of a new boson at the LHC mi, whose properties are compatible with those 
of the SM Higgs. They are also important in many new physics scenarios 13 [IQ], in particular in 
the search for extra-dimensions m or supersymmetry [ 12 ]. 

In this paper we are interested in the process pp —?■ 77 W, and in particular in the transverse- 
momentum [qt) spectrum of the diphoton pair. The lowest-order process ( 0 ( 05 )) occurs via the 
quark annihilation subprocess qq —)■ 77 . The QCD corrections at the hrst order in the strong 
coupling as involve quark annihilation and a new partonic channel, via the subprocess qg —)■ 77 ^. 
First order corrections have been computed and implemented in several fully-differential Monte 
Carlo codes HQi Emsi do]. At the second order in the strong coupling as the gg channel starts 
to contribute, and the large gluon-gluon luminosity makes this channel potentially sizeable. 

The amplitudes needed to evaluate the corrections at the second order in the strong coupling 
Q! 5 , for diphoton production, have been presented in (TT] [THl dQ] ED], and hrst put together in 
a complete and consistent 0{a'^) calculation in the 27 NNLD code [21]. The next-order gluonic 
corrections to the box contribution (which are part of the N^LO QCD corrections to diphoton 
production) were also computed in ref. [H] and found to have a moderate quantitative effect. 

The calculation of the qx spectrum poses an additional challenge with respect to more inclu¬ 
sive calculations, such as the total cross section. In the large-g^ region ~ M^^), where the 
transverse momentum is of the order of the diphoton invariant mass calculations based on 
the truncation of the perturbative series at a hxed order in as are theoretically justihed. In this 
region, the QCD radiative corrections are known up to the next-to-leading order (NLO), including 
the corresponding partonic scattering amplitudes with X = 2 partons (at the tree level [18]) and 
the partonic scattering amplitudes with X = 1 parton (up to the one-loop level [ID])- We remind 
the reader that at least one additional parton is needed in order to have gr 7 ^ 0 for the diphoton 
pair. The qx spectrun of the diphoton pair has been calculated in fully-differential Monte Carlo 
codes at LO [IS] EU ES] ED] and at NLO [HI ED] ED]- Recently, hrst calculations for diphoton 
production in association with two [DU ESI ED] and three [26] jets at NLO became available. 

The bulk of the diphoton events is produced in the small-gT region (g^ <C M^^), where the 
convergence of the hxed-order expansion is spoiled by the presence of large logarithmic terms, 
ag ln"*(M^^/g|,). In order to obtain reliable predictions these logarithmically-enhanced terms have 
to be systematically resummed to all perturbative orders [2^-[38]. The resummed calculation, 
valid at small values of g-r, and the hxed-order one at large qx have then to be consistently matched 
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to obtain a pQCD prediction for the entire range of transverse momenta. 

We use the transverse-momentum resummation formalism proposed in Refs. [3711391 HD] (see 
also [H] for processes initiated by gg annihilation). The formalism is valid for a generic process 
in which a high-mass system of non strongly-interacting particles is produced in hadron-hadron 
collisions. The method has so far been applied to the production of the Standard Model (SM) 
Higgs boson [391 soma Hal E], Higgs boson production in bottom quark annihilation [15] , Higgs 
boson production via gluon fusion in the MSSM [16], single vector bosons at NLL-I-LO m and 
at NNLL-I-NLO [H], WW [IHl |50] and ZZ [51] pairs, slepton pairs [52], and DY lepton pairs in 
polarized collisions [53] - 

Finally, note that besides the direct photon production from the hard subprocess, photons can 
also arise from the fragmentation of QCD partons. The computation of fragmentation subprocesses 
requires (the poorly known) non-perturbative information, in the form of parton fragmentation 
functions of the photon (the complete single- and double-fragmentation contributions are imple¬ 
mented in DIPHOX [13] for diphoton production at the hrst order in as). However, the effect of 
the fragmentation contributions is sizeably reduced by the photon isolation criteria that are nec¬ 
essarily applied in hadron collider experiments to suppress the very large irreducible background 
{e.g., photons that are faked by jets or produced by hadron decays). Two such criteria are the so- 
called “standard” cone isolation and the “smooth” cone isolation proposed by Frixione [5l]- The 
standard cone isolation is easily implemented in experiments, but it only suppresses a fraction 
of the fragmentation contribution. By contrast, the smooth cone isolation (formally) eliminates 
the entire fragmentation contribution. For all of the results presented in this paper we rely on 
the smooth isolation prescription, which, for the parameters used in the experimental analysis 
reproduces the standard result within a 1% accuracy [55] . 

The paper is organized as follows. In Sect. |2]we briefly review the resummation formalism of 
Refs. [371 [39l 00]. In Sect. [3] we present numerical results and we comment on their comparison 
with the LHC data [6]. We also study the scale dependence of our results with the purpose of 
estimating the corresponding perturbative uncertainty. In Sect. 0]we summarize our results. 


2 Transverse-momentum resummation 


We briefly recall the main points of the transverse-momentum resummation formalism of Refs. m 
l39| 00], referring to the original papers for the full details. The formalism is general, as long as 
the measured hnal state is composed of non strongly-interacting particles (transverse-momentum 
resummation for strongly-interacting hnal states, such as heavy-quark production, has been de¬ 
veloped in Refs. [561157]). Here we specialize to the case of diphoton production only for ease of 
reading. The inclusive hard-scattering process considered is 

hi{pi) + h2{p2) + X , (1) 

where hi and h 2 are the colliding hadrons with momenta pi and p 2 , 77 is the diphoton pair 
with invariant mass transverse momentum and rapidity ?/, and X is an arbitrary and 

undetected hnal state. 

The corresponding fully diherential cross section, in q^, and y, which we denote for simplic- 
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ity (since our focus is on the qt distribution) by da-y-y/dq^, can be written using the factorization 
formula as 


da, 


77 


dq^ 




Myy, S) 


a,b 


^^77 

1/^2 fa /hi (^1 )/^f) /v^2 (^2) ^^2 idTi ^'y'yi y 1 S, Oig, 


' yV) 


( 2 ) 


(up to power-suppressed corrections), where the fa/h{x,jfj^) (a = q,q,g) are the parton densities 
of the hadron h at the factorization scale /if, as = daff/dqf is the pQCD partonic cross 

section, s {s = xiX 2 s) is the square of the hadronic (partonic) centre-of-mass energy, and /i^ is 
the renormalization scale. 


In the region where qr ~ Myy the QCD perturbative series is controlled by a small expansion 
parameter, as{Myy), and a hxed-order calculation of the partonic cross section is theoretically 
justihed. In this region, the QCD radiative corrections are known up to next-to-leading order 

(NLO) [niiiHiiisiED]. 

In the small-gT region (g^ "C Myy), the convergence of the hxed-order perturbative expansion 
is spoiled by the presence of powers of large logarithmic terms, ag ln™'(M.^^/g|.). To obtain reliable 
predictions these terms have to be resummed to all orders. 


To perform the resummation, we start by decomposing the partonic cross section as 


d»i2 

dq^ dq^ dq^ 


(3) 


The hrst term on the right-hand side contains all the logarithmically-enhanced contributions, 
which have to be resummed to all orders in as, while the second term is free of such contributions 
and can thus be evaluated at hxed order in perturbation theory. Using the Fourier transformation 
between the conjugate variables qx and b {h is the impact parameter), the resummed component 
da^^^^2 can be expressed as 


^^■yyab , 2 2\_ ^^“^77 

^^2 \QTiMyy,y,S,0!g,p,x,,p,p) ^ 


db - Joibgr) {b, Myy, y, s; as, ^r, yl) , (4) 


where Jo{x) is the Oth-order Bessel function. The form factor is best expressed in terms of 
its double Mellin moments taken with respect to the variables Zi, at hxed Myy, with 


Z1Z2 = 2: = 



^2 


the resummation structure of can be organized in an exponential form [j 


(5) 


y^JhN2(b, Myy, y, as, ^r, fip) — {Myy, as; M^^/y,p, M^^/fip, M^^/ 

xexp{gN,N2ias,L;M^,^/y%M^,^/l4^^)} , ( 6 ) 

were we have dehned the logarithmic expansion parameter L = In^yl^gb"^/b^), and bo = {■jp = 

0.5772... is the Euler number). The scale yres {yres ~ Myy), which appears on the right-hand side 

^For the sake of simplicity we consider here only the case of the diagonal terms in the flavour space of the 
partonic indices a, b. For a detailed discussion, we refer to Ref. [SillS]- 
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of Eq. ([2]), is the resummation scale [3S]. Although (i.e., the product x expl^ATjATj}) 

does not depend on when evaluated to all perturbative orders, its explicit dependence on 
fires appears when is computed by truncation of the resummed expression at some level 

of logarithmic accuracy (see Eq. ([7]) below). Variations of fires around can thus be used to 
estimate the size of yet uncalculated higher-order logarithmic contributions. 

The form factor expl^TVjATj} is universalil and contains all the terms that order-by-order in as 
are logarithmically divergent as 6 —)■ cxo (or, equivalently, qT —)■ 0). The resummed logarithmic 
expansion of the exponent GniN 2 is dehned as follows: 

^Vi7V2(as, L; fij^, fiU = L g^^\asL) + M^^/fi% M^^/fiU 

+ ~fi'WiV2(“S-^, ^77/hL ^^'r/^^les) + • • • (’^) 

where the term L g^^^ collects the leading logarithmic (LL) contributions, the func- 

tion fi']viV 2 includes the next-to-leading leading logarithmic (NLL) contributions [32], 

9nIn2 controls the NNLL C>(af+"L"“^) terms [33l[3ll[58l|59| and so forth; p is the number of pow¬ 
ers of as in the LO (Born) process. In Eq. (f^, a^L is formally of order 1, so there is an explicit 
0{as) suppression between different logarithmic orders. The explicit form of the functions g^^\ 
9 nIn2 9nIn2 found in Ref. [39] . The process dependent function does not depend 

on the impact parameter b and it includes all the perturbative terms that behave as constants as 
6 —)■ cxo. It can thus be expanded in powers of as: 






[^y^VSHKIA-, KIA, Ml/p-l.) 


where is the partonic cross section at the Born level. Since the formalism applies to non 
strongly-interacting hnal states, in general the Born cross-section can only correspond to a qq or 
gg initial state. In the specihc case of the diphoton production, both channels contribute, but at 
different orders in a*: the qq subprocess initiates as a pure QED process (C>(as)°), while the gg 
one requires a fermion loop, starting at 0{as)‘^- 


In the present work, we keep contributions up to an uniform order in a^ (and all orders in 
asT), namely up to agL"'~^. For the qq channel, this requires the inclusion of the "H coefficients 
of Eq. ([8]) up to order 2: the hrst-order coefficients 'Hjf^N 2 known since a long time [58], while 
the second-order coefficients were computed only relatively recently [7T], |3B] • For the gg 

channel, it is sufficient to include the leading l-i contribution (that is, the Born cross-section) and 
the appropriate Q in the exponential of Eq. d?]). Since it does not require any additional numerical 
effort, we decided, in all the plots presented in the paper, to include all the terms up to in the 
exponential Q factor also for this channel. In this way, we technically include some terms which 
are of higher order in a^ with respect to those in the qq channel; however we checked that those 
terms result in a negligible numerical effect (at 1% accuracy), that is, the difference produced by 
including the higher order terms is within the error bands obtained by the scale variations, which 
verihes the stability of the calculation. 

^The form factor does not depend on the final state; all the hard-scattering processes initiated by qq {gg) 
annihilation have the same form factor. 
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Within a straightforward (‘naive’) implementation of Eq. (E]), the resummation of the large 
logarithmic contributions would affect not only the small-gT region, but also the region of large 
values of qt- This can easily be understood by observing that the logarithmic expansion parameter 
L diverges also when & —)■ 0. To reduce the impact of unjustihed higher-order contributions in the 
large-gr region, the logarithmic variable L in Eq. ([6]) is actually replaced by L = In -|- 1) 

[391112]. This (unitarity related) replacement has an additional and relevant consequence: after 
inclusion of the hnite component (see Eq. ([9])), we exactly recover the hxed-order perturbative 
value of the total cross section upon integration of the qt distribution over qt (i.e., the resummed 
terms give a vanishing contribution upon integration over qt)- 

We now turn to consider the hnite component of the transverse-momentum cross section (see 
Eq. (El)). Since does not contain large logarithmic terms in the small-gT region, it can be 

evaluated by truncation of the perturbative series at a given hxed order. In practice, the hnite 
component is computed starting from the usual hxed-order perturbative truncation of the partonic 
cross section and subtracting the expansion of the resummed part at the same perturbative order. 
Introducing the subscript f.o. to denote the perturbative truncation of the various terms, we have: 


7 ^ (fin.) 




7 ^ (res.) 

- dqf . 

f.o. 

. dqf . 

f.o. 

. dqf . 


( 9 ) 


This matching procedure between resummed and hnite contributions guarantees to achieve uni¬ 
form theoretical accuracy over the region from small to intermediate values of transverse momenta. 
At large values of qt, the resummation (and matching) procedure is eventually superseded by the 
customary hxed-order calculations (their theoretical accuracy in the large-gr region cannot be 
improved by resummation of the logarithmic terms that dominate in the small-g^ region). 

In summary, the inclusion of the functions 9nIn2^ resummed component, 

together with the evaluation of the hnite component at LO (i.e. at 0(as)), allows us to perform 
the resummation at NLL-I-LO accuracy. This is the theoretical accuracy used in previous studies 
[T6| IHn] [6T| [62] of the diphoton qt distribution. Including also the functions and 

together with the hnite component at NLO (i.e. at 0{ag)) leads to full NNLL-I-NLO accuracy. 

Using the coefficient [2H 138] . we are thus able to present the complete result for 

the diphoton g^-distribution up to NNLL-I-NLO accuracy. We point out that the NNLL-I-NLO 
(NLL-I-LO) result includes the full NNLO (NLO) perturbative contribution in the small-gy re¬ 
gion. In particular, the NNLO (NLO) result for the total cross section is exactly recovered upon 
integration over qT of the diherential cross section da^^/dqx at NNLL-I-NLO (NLL-I-LO) accuracy. 


We conclude this section with some comments on the numerical implementation of our calcu¬ 
lation. Within our formalism, the resummation factor is directly dehned, at hxed 

in the space of the conjugate variables b and Ni, N 2 . To obtain the hadronic cross section, 
we have to perform inverse integral transformations: the Bessel transformation in Eq. ((1]) and 
an inverse Mellin transformation, implemented following the prescription introduced in Ref. [63] . 
These integrals are carried out numerically. The Mellin inversion requires the numerical evalu¬ 
ation of some basic Wmoment functions that appear in the expression of the the second-order 
coefficients [2ll 13H]: this evaluation has to be performed for complex values of Ni and W; 

to evaluate some of the needed special function at complex value, we use the numerical routines 
of Ref. [M]- We recall [39] that the resummed form factor , L)} is singular at 
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the values of b where q;s(/U^)L > 7r//9o {/3q is the hrst-order coefficient of the QCD (3 function). 
We avoid this singularity by introducing a smooth effective cut-off at small b, which is shown to 
have a negligible effect in the final result. 

It is known that at small values of qt, the perturbative QCD approach has to be supplemented 
with non-perturbative contributions, since they become relevant as qt decreases. A discussion on 
non-perturbative effects on the qr distribution is presented in Ref. [39], and related quantitative 
results are shown in Sect. [3l 


3 Numerical results for diphoton production at the LHC 


In this section we consider diphoton production in pp collisions at LHC energies = 7 TeV). 
We present our resummed results at NNLL-I-NLO accuracy, and compare them with NLL-I-LO 
predictions and with available LHC data [B]. Since the present formulation of the qx resummation 
formalism, is restricted to the production of colourless systems F, it does not treat parton frag¬ 
mentation subprocesses (here F includes one or two coloured partons that fragment); therefore, 
we concentrate on the direct production of diphotons, and we rely on the smooth cone isolation 
criterion proposed by Frixione la (see also Ref. [65] [66]) which is defined by requesting 

^ < Ft max x(r) , 

inside any r'^ = {y — + (0 ~ > (10) 


with a suitable choice for the function x(r). This function has to vanish smoothly when its 
argument goes to zero 0 , if r ^ 0), and it has to verify 0 < x(r) < 1, if 0 < r < i? . 

One possible choice is 

/I — cos(r)\ "■ 

\ 1 — cos R ) ’ 

where n is typically chosen as n = 1. This condition implies that, closer to the photon, less 
hadronic activity is allowed inside the cone. At r = 0, when the parton and the photon are 
exactly collinear, the energy deposited inside the cone is required to be exactly equal to zero, 
and the fragmentation component (which is a purely collinear phenomenon in perturbative QCD) 
vanishes completely. Since no region of the phase space is forbidden, the cancellation of soft gluon 
effects takes place as in ordinary infrared-safe cross sections. That is the main advantage of this 
criterion: it eliminates all the fragmentation component in an infrared-safe way. By contrast, it can 
not be implemented within the usual experimental conditions; the standard way of implementing 
isolation in experiments is to use the prescription of Eq. fITT]) with a constant x(^) = 1- In any 
case, from a purely pragmatic point of view, it has been recently shown [55] that if the isolation 
parameters are tight enough (e.g.. Ft max < 6 GeV, R = 0.4), the standard and the smooth cone 
isolation prescription coincide at the 1% level, which is well within the theoretical uncertainty of 
our predictions. 



The acceptance criteria used in this analysis {y/s = 7 TeV) are those implemented by the 
ATLAS collaboration analysis [6]; in all the numerical results presented in this paper, we require 
^harder > 25 GeV, > 22 GeV, and we restrict the rapidity of both photons to satisfy 

\y^\ < 1.37 and 1.52 < \y^\ < 2.37. The isolation parameters are set to the values Ft max = 4 GeV, 
n = 1 and R = 0.4, and the minimum angular separation between the two photons is R^^ = 0.4. 
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We use the Martin-Stirling-Thorne-Watt (MSTW) 2008 [HZ] sets of parton distributions, with 
densities and as evaluated at each corresponding order (i.e., we use (? 7 , + l)-loop as at N"LO, with 
n = 0,1,2), and we consider Nf = 5 massless quarks/antiquarks and gluons in the initial state. 
The default renormalization (/i/j) and factorization (fip) scales are set to the value of the invariant 
mass of the diphoton system, Hr = Hf = M 77 , while the default resummation scale (pres) is set 
to Hres = M^^/2. The QED coupling constant a is hxed to a = 1/137. 




Figure 1: The qt spectrum of the photon pair (solid red lines) at the LHC (1 TeV): results 
at NLL+LO (left panel) and NNLL+NLO (right panel) accuracy. Each result is compared to 
the corresponding fixed-order result (dashed lines) and to the finite component (dotted lines) in 
Eg. 1^. The resummed spectrum includes a non-perturbative (NP) contribution parametrized as 
in Eg. oa and is obtained within the 3-body approach. 


This choice of the order of the parton densities and as is fully justihed both in the small-g^ 
region (where the calculation of the partonic cross section includes the complete NNLO (NLO) 
result and is controlled by NNLL (NLL) resummation) and in the intermediate-gr region (where 
the calculation is constrained by the value of the NNLO (NLO) total cross section). 

Non-perturbative (NP) effects are expected to be important at very small qt- In this paper we 
follow the strategy of Ref. [39], implementing them by multiplying the 6 -space form factor 
of Eq. dH) by a ‘NP factor’ which consists of a gaussian smearing of the form 

S%p = expi-CagNph"^), ( 12 ) 

where a denotes the initial state channel, a = F for qq and a = A for gg (as usual, Cp = 
(A(? — l)/(2iVc) and Ca = W)- In order to asses the importance of the NP contributions, we vary 
gNP in the interval from g^p = 0 GeV^ (no NP contributions) to gNP = 2 GeV^, corresponding 
to moderate NP effects [39] . 

An additional and potentially important source of theoretical uncertainty arises from an ambi¬ 
guity in the dehnition of the photon momenta in the resummation formalism. In fact, in the main 
resummation formula (jl]), which is used to dehne both the resummed and (via the subtraction 
Eq. ([9|)) the hnite contributions to the partonic cross-section, the diphoton pair total transverse 
momentum qp is not associated with the recoil of any extra physical particles in the hnal state. 
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After qx resummation the angular distributions of the photons are still provided by the Born 
level functions which appear as multiplicative factors in front of the Sudakov form 

factor of Eq. ([7]). At this point there are two strategies to follow, which differ by corrections that 
are of 0{qT/ order-by-order in the perturbative expansion j69] (after having matched the 
resummed calculation with a complete N^LO calculation, these corrections start to contribute at 
the N^+^LO). 




Figure 2: Uncertainties in the qr spectrum of diphoton production at the LHC (y/s = 7 TeV). In the 
left panel we show the different predictions associated with the 2-hody and 3-hody parametrization 
of the photon momenta at the central scale pp = Ur = p^es = M^^/2. In the right panel, 
the bands (solid and dotted lines) are obtained by varying qip and fiR as described in the text; the 
prediction at the central scale with the NP contribution turned off (g^p = 0 GeV^ in Eq. ca; is 
also shown. 

One of them is to just use the Born phase space for all the angular distributions in front of 
the Sudakov form factor of Eq. ([7]). In this case the transverse momentum qp is neglected in the 
calculation of the photon momenta, while the momenta g,, {i = 1, 2) of the colliding partons are 
given by: 

qf = XiPf qlf = X 2 Pf , (13) 

where Xi {i = 1, 2) are the parton momentum fractions and Pi {i = 1, 2) are the momenta of the 
colliding hadrons. The momenta g* respects the Born level kinematics 

gi + g 2 = q^i + g 72 , (14) 

where q^i , (z = 1, 2) are the momenta of the two photons in the hnal state. In this case the two 
photons are always in the back-to-back conhguration, and therefore the kinematic effects of the 
transverse-momentum recoil, are not included in the momentum of each single photon. We call 
this approach the 2-body phase space; all the differential distributions which use this method have 
a (2-body) label. 

A more elaborate, and arguably more physical, approach is to consider the effects of the 
transverse-momentum recoil in the two-photon hnal state. Therefore these kinematic effects have 
to be ‘absorbed’ by the incoming parton momenta gi and g 2 , in order to respect the Born level 
kinematics 


9 i ( 9 t ) + q^iqr) = + q^iiqr) , 


(15) 






and the LO kinematics of Eq. fll3p . There are different consistent implementations of this ap¬ 
proach, all of them giving equivalent (up to higher order corrections) results [§. The implementa¬ 
tion that we use in the following is essentially equivalent to dehne the separate photon momenta 
according to the Born level angular distribution computed in the Collins-Soper frame isn of the 
diphoton system (in practice, this is the same procedure used in Ref. USD- We call this method 
the 3-body phase space, and distributions obtained using this method are labelled (3-body). The 
modihcations introduced by this phase space parametrization in the momentum fractions Xi, X 2 
of the incoming partons are neglected because they produce negligible effects {0{qT/M^^)) in the 
cross-section. Notice that it is not the formal 3-body phase space (the 3-body effect is produced by 
the Lorentz transformation with hnite qt of the Collins-Soper frame). We have checked that if the 
formal 3-body phase space is used0, the changes in the cross section (comparing the 2-body and 
3-body cases) are of order 0{qT/M^^). The ambiguity in the treatment of the photon momenta 
is considered as an additional source of theoretical uncertainty. 



0 5 10 15 20 25 30 35 40 


qTvv [GeV] 



Figure 3: The qx spectrum of diphoton production at the LHC (^/s = 7 TeV). Here, the scales 
Tf = Tr = o,re kept fixed while we vary the resummation scale pires to establish its impact on 
the cross section. In the left panel we show the range 0 GeV< ^^77 < 40 GeV, and in the right 
panel the full spectra in logarithmic scale. 


In Fig. m left panel, we present the NLL-I-LO qx spectrum at the LHC (i/s = 7 TeV). The 
NLL-I-LO result (solid line) at the default scales (pf = hJ? = M 77 ; Pres = ^ 77 / 2 ) is compared with 
the corresponding LO result (dashed line). We use the 3-body phase space and a NP parameter 
gMP = 2 GeV^. The LO hnite component of the spectrum (see Eq. ([3])), is also shown for 
comparison (dotted line). We observe that the LO result diverges to -|-cxo as qx —)■ 0, as expected. 
The hnite component is regular over the full qx range, it smoothly vanishes as t 0 and 
gives an important contribution to the NLL-I-LO result in the low-g^ region. That is mostly 
originated by the qg channel, which starts at NLO and provides a subleading correction in terms 
of logs (single logarithmic terms) but contributes considerably to the cross-section due to the 
huge partonic luminosity compared to the formally leading qq channel. The resummation of the 

^The details of the implementation of these kinematic effects are discussed in a forthcoming paper [69) . 

^Which violates the Born level kinematics of Eq. (ITTl) . 
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small-gr logarithms leads to a well-behaved distribution; it vanishes as ^ 0 and approaches 
the corresponding LO result at large values of qt- 

The results in the right panel of Fig. [1] are systematically at one order higher: the qt spectrum 
at NNLL-I-NLO accuracy (solid line) is compared with the NLO result (dashed line) and with the 
NLO hnite component of the spectrum (dotted line). The NLO result diverges to —cxo as qt ^ 0 
and, at small values of qt, it has an unphysical peak that is produced by the compensation of 
negative leading and positive subleading logarithmic contributions. The contribution of the NLO 
hnite component to the NNLL-I-NLO result is of the order of the 50% at the peak and becomes 
more important as qt increases. A similar quantitative behaviour is observed by considering 
the contribution of the NLO hnite component to the NLO result. At large values of the 
contribution of the NLO hnite component tends to the NLO result. This behaviour indicates that 
the logarithmic terms are no longer dominant and that the resummed calculation cannot improve 
upon the predictivity of the hxed-order expansion. 

We also observe that the position of the peak in the NNLL-I-NLO qt distribution is slightly 
harder than the corresponding NLL-I-LO qt distribution. This effect is (in part) due to the large 
transverse-momentum dependence of the hxed order corrections. 

As discussed in Sect. [H the resnmmed calculation depends on the factorization and renormal¬ 
ization scales and on the resnmmation scale fires- Onr convention to compute factorization and 
renormalization scale uncertainties is to consider independent variations oi fip and fiR by a factor 
of two around the central values fip = fJ-R = Af .77 in independent way in order to maximise them: 
[fip = 2 My.y, fiR = M^^/2, fires = ^ 77 / 2 ) and {fiR = 2 fip = M^^/2, fires = M^-y/2). The 
uncertainty due to the resummation scale variation is assessed separately by varying it between 
fires = and fires = M^-y at hxed fip and fiR. 



Figure 4: The qt spectrum of diphoton pairs at the LHC. The NNLL+NLO result is compared with 
the NLL+LO result, for the window 0 GeV < g^ < 40 GeV (left panel) and the full spectra (right 
panel). We use the 3-body parametrization in the resummed cross-section, and set = 2 GeW. 
The bands are obtained by varying fiR and fip as explained in the text. 

In Fig. [2] we compare the impact of the various sources of theoretical error for the NNLL-I-NLO 
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Figure 5: Comparison of the theoretical prediction for the qt spectrum of diphoton pairs at the 
LHC to the experimental data. The NNLL+NLO result is compared with the ATLAS data of 
Ref JBS, for the window 0 CeV < gr < 40 CeV (left panel) and the full spectra (right panel). 
In the right panel the NNL+LO distribution at central scale (dotted line) is also shown in order 
to compare it with the data and the NNLL+NLO result. In the theoretical curves we use the 2- 
body parametrization for the resummed cross-section and set g^p = 2 Ce\R; the bands (solid and 
dashed lines) are obtained by varying pr and as explained in the text. 

predictions: the 2-body and 3-body parametrizations of the photon phase space, the effect of 
the variation of the non-perturbative contribution, and the variation of the factorization and 
renormalization scales. In the left panel of Fig. [2] we use only the central scale {pp = Pr = 

Pres = M.y.y/2), and we note that the 3-body method results in a slightly larger cross section 
(about 10%) around the peak {qp ~ 5 GeV) than the 2-body one. For qp > 20 GeV all the 
contributions in the left panel coincide, which is consistent with the fact that at these values the 
resummed component starts to vanishes and the hxed order result dominates the cross section, 
as we can anticipate from Fig. [T] Also we note that if we use a NP parameter, the peak of the 
distribution is located at larger values oi qp {qp ^ 5 GeV) and the shape of the distribution is 
slightly different (for values 0 GeV< qp < 22 GeV) from the case in which the NP parameter is 
not implemented. These differences, which are stronger for qp < 10 GeV, have their origin in the 
resummed component, which is the only contribution that depends on gpfp- In the right panel of 
Fig. [2] we show the comparison between the variation of the scaled in the 3-body approach with 
the central scale results of the 2-body and 3-body frameworks (with and without non perturbative 
parameter, respectively). Evidently, the uncertainty due to the ambiguity of the parametrization 
of the photon momenta turns out to be subdominant with respect to the one arising from the 
variation of the scales which provides, by itself, a reasonable estimate of theoretical uncertainties. 

In Fig. [3] we show the NNLL-I-NLO transverse momentum distribution for three different 
implementations of the pres parameter {pres = at hxed Pp = Pr = M.y.y. 

The 2-body phase space and a non-perturbative parameter g^p = 2 GeV^ were used. We notice 
the small impact of the variation of the pres scale in the cross section (at per-cent level). In the 
left panel of Fig. [3] we present the transverse momentum distribution for values of qp within the 

II We vary the scales from [yp = 2M.y^\yR = Mj.y/2) to {rr = 2Mj.y-, yp = Mry.y(2) at fixed yres = M.y..^/2. 
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Figure 6: Comparison of the theoretical prediction for the qt spectrum of diphoton pairs at the 
LHC with the experimental data. The NNLL+NLO result is compared with the ATLAS data of 
Ref. 0, for the window 0 CeV < gr < 40 CeV (left panel) and the full spectra (right panel). In 
the right panel the NNL+LO distribution at central scale (dotted line) is also shown in order to 
compare it with the data and the NNLL+NLO result. In the theoretical curves we use the 3-body 
parametrization for the resummed cross-section and set g^p = 2 CeV^. The bands (solid and 
dashed lines) are obtained by varying /i^ and as explained in the text. 

interval 0 GeV< qp < 40 GeV, and in the right panel of Fig. |3]the full spectra. We also notice 
that the strongest effect of the variation of the pres scale appears in the last bin of right panel 
of Fig. [3l This is expected since the resummation scale effectively sets the value of transverse 
momentum at which the logarithms are dominant. A choice of a very large resummation scale 
affects the distribution at larger transverse momentum and might in general result in a mismatch 
with the fixed order prediction due to the artificial introduction of unphysically large logarithmic 
contributions in that region. Similar results are obtained if the 3-body phase space is used instead 
of the 2-body one. 

In Fig. m we compare the variation of the scales of the NNLL-I-NLO and NLL-I-LO predictions 
(3-body phase space), for the interval 0 GeV< qp < 40 GeV (left panel) and the full spectra 
(right panel). We notice that the dependence on the scales is not reduced when going from 
NLL-I-LO to NNLL-I-NLO. This is mostly because at NNLL-I-NLO a new channel (gg) opens, in 
which the box contribution (effectively “LO” but formally (9(a|)) ruins the reduction of the scale 
dependence usually expected when adding second order corrections for the qq channel and first 
order corrections for the qg channel. This effect has the same origin that the reported behaviour 
of the diphoton production at NNLO of Ref. [21], when the variation of the scales is implemented. 
Since NNLL-I-NLO is the hrst order at which all partonic channels contribute, it is possible to 
argue that this is the hrst order at which estimates of theoretical uncertainties through scale 
variations can be considered as reliable. The same results are obtained if the 2-body phase space 
is used instead of the 3-body one. 

The peak observed in the right panel of Fig. IHis the so called Cuillet shoulder [6H], which is a 
real radiation effect and has its origin in the hxed order contribution. It appears stronger in the 
NNLL-I-NLO qp distribution than in the NLL-I-LO, due to the larger size of the real contributions 
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Figure 7: Individual photon spectra in diphoton production at the LHC. Left panel: the spec¬ 
trum of the hard photon at the LHC JEi (y/s = 7 TeV). Right panel: the spectrum of the softer 
photon at the LHC (y/s = 7 TeV). In both panels, we compare the resummed prediction using 
the 3-body parametrization in the resummed cross-section (gNP = 2 Ce\R) with the fixed order 
prediction. 

at NLO. 

In the small gr region (g^ < 4 GeV) the real radiation effects are no longer dominant in the 
gr distribution. The finite and resummed component (which vanish as qt —t 0) are of the same 
order for gj^ < 4 GeV. In absence of the strong real radiation effects the contributions are almost 
completely Born like. This is the main reason why the NNLL+NLO bands overlaps with those at 
the previous order. 

In Figs. Eland[H] we compare the LHG data (i/s = 7 TeV) from ATLAS [6] with our resummed 
theoretical predictions (at NNLL+NLO and NLL+LO) using the 2-body and 3-body approaches, 
respectively. In both cases we use a NP parameter different from zero (gwp = 2 GeV^) and we 
estimate the theoretical uncertainty by the variation of the and fip scales. In the left panels 
we show the g^ distribution in the window (0 GeV< g^ < 40 GeV), while in the right ones we 
show the full spectra in logarithmic scale. 

We observe in general an excellent agreement between the resnmmed NNLL+NLO prediction 
and the experimental data, that is accnrately described within the theoretical uncertainty bands 
in the whole kinematic range. Also we observe that the NLL+LO result is not enough to describe 
the phenomenology of the transverse-momentum distribution of the LHG data (Figs. [5] and [6l right 
panel). By direct comparison to the hxed order prediction, we notice that the effect of resummation 
is not only to recover the predictivity of the calculation at small transverse momentum, but also 
to improve substantially the agreement with LHG data [H]. 

While the resnmmation performed in this work reaches NNLL accnracy formally only for the 
diphoton transverse momentum distribution, its predictions can be extended to other observables 
as well, since at least the leading logarithmic contributions have a common origin from soft and 
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collinear emission. Note that the 3-body approach is more suitable for a consistent implementation 
of these leading logarithmic effects on the observables. In Fig. [3 we show results on more exclusive 
observables: the qt distributions of the harder (left-hand plot) and softer (right-hand plot) photon 
at the central scale. We compare the results obtained with the 2-body and 3-body phase space 
approaches, using a non-perturbative parameter g^p = 2 GeV^ with the hxed order result at 
NNLO for diphoton production 121 ]. 

The 2-body phase space transverse-momentum distribution at NNLL-I-NLO provides the same 
result than the NNLO fixed order cross section for diphoton production |2T]. This is consistent 
with two following related facts: i) the single photon momentum does not carry any information 
about the recoil due to the transverse momentum in the 2-body approach; ii) the NNLO (NLO) 
result for the total cross section is exactly recovered upon integration over qt of the differential 
cross section da.^^/dqT at NNLL-I-NLO (NLL-I-LO) accuracy. 

The effects of resummation are only present in these more exclusive observables if the recoil due 
to qT are absorbed by the photons in the final state 0, which is equivalent to the implementation 
of a 3-body like phase space. In this way, because the single photon momentum depends on q^p, the 
integral over qx of the NNLL-I-NLO (NLL-I-LO) distribution does not recover the NNLO (NLO) 
result as we can observe in Fig. [3 Here the last integral over the single photon momentum of the 
NNLL-I-NLO (NLL-I-LO) distribution, is required in order to recover the NNLO (NLO) result for 
the total cross section. 

We comment on the qp distribution of the softer photon in the region around the back-to-back 
threshold q^^ ~ 25 GeV (see Fig. [3 right panel). The NLO fixed order result has a step-like be¬ 
haviour, and this necessarily produces ITU integrable logarithmic singularities at each subsequent 
perturbative order. The peak of the NLO hxed order distribution at q^ ^ ^ 25 GeV is an artifact 
of these perturbative instabilities. The instability is cured by all-order perturbative resummation, 
which leads to a smooth qp distribution with a shoulder-like behaviour lu in the vicinity of the 
back-to-back threshold. 

In Fig. [8] we present the results of the cross section as a function of the azimuthal angle 
In the left panel we compare the hxed order (NLO), hnite (NLO) and full (NNLL-I-NLO) 
distributions. The hxed order component dominates the cross section over the whole A^..^..^ range. 
However, as could be expected, the ehect of resummation is stronger for kinematic conhgurations 
near the A<h..y..y ~ vr which correspond to gr ~ 0 GeV. As in the case of the hxed order qp 
distribution, the A<h..y..y hxed order diherential cross section is not well-behaved near the back-to- 
back conhguration: it actually diverges as A<F.^.^ vr {qp 0). The hnite contribution (Eq. (jU])) 
is well-behaved near the back-to-back conhguration, and the full result (NNLL-I-NLO) improves 
the description in the region near A$.y.y ~ 0. 

In the right panel of Fig. [8] we compare our theoretical prediction at NNLL-I-NLO level of 
accuracy with the LHG data [6] using the variation of the fip and /if scales to estimate the 
theoretical uncertainty. We observe that the transverse momentum resummation provides a better 
description of the data with respect to the hxed order result. In both panels of Fig. [H] we used 
the 3-body approach to describe the diphoton phase-space. In fact, in this kind of observables 
(also see Fig. (3), the 2-body parametrization again reproduces the result of the hxed order cross 

**And also the recoil due to qp has to be absorbed in the initial state to restore the momentum conservation 
(see Eq. (fTSI l. 
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Figure 8: The A$-y-y distribution of diphoton pairs at the LHC. In the left panel, we show the 
fixed order prediction, the complete resummed prediction and just the resummed contribution at 
the central scale. In the right, the full NNLL+NLO result (using the extreme values for and 
Pp to estimate the theoretical uncertainty) is compared with the ATLAS data of Ref. 


section. 


4 Summary 


In this paper we performed the transverse momentum resummation for diphoton production at 
NNLL accuracy in hadron collisions. At small values of qp, the calculation includes the resumma¬ 
tion of all logarithmically-enhanced perturbative QCD contributions, up to next-to-next-to-leading 
logarithmic accuracy; at intermediate and large values of qp, it combines the resummation with 
the fixed next-to-leading order perturbative result. The combination is performed in such a way 
as to exactly reproduce the known next-to-next-to-leading order result for the total cross section; 
in the end, the calculation consistently includes all perturbative terms up to formal order The 
theoretical uncertainty was estimated by varying the various scales (renormalization, factorization 
and resummation) introduced by the formalism as well as the parametrization of the diphoton 
phase-space. The result was compared to experimental data, showing good agreement between 
theory and experiment over the whole qp range. With respect to the fixed-order calculation, the 
present implementation provides a better description of the data and recovers the correct physical 
behaviour in the small qp region, with the spectrum smoothly going to zero. The same set-up also 
allows the calculation of more exclusive observable distributions; the qp spectrum of the individual 
photons and the distribution are given as examples. 
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